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Abstract 

TD ■ 

<—j ■ In this paper, we construct a class of spin-1/2 antiferromagnetic (AFM) 

Q ' two-chain ladder models consisting of blocks of n-spin tetrahedral clusters 

O , alternating with two-spin rungs. For n=4 and 6 and in extended parame- 
ter regimes, the exact ground state of the ladder is shown to be a product 

fvj ' of the ground states of the rungs and the n-spin blocks, in both zero and 

^. , finite magnetic fields. In the latter case, magnetization/site m versus ma- 

£C) ' genetic field h plot exhibits well-defined magnetization plateaus. 

Q\ ' PACS numbers: 75.10 Jm, 75.40 Mg, 75.50Ee 

m ! 

r — , Spin ladders have been widely studied in recent times as these systems exhibit 

a variety of novel phenomena in the undoped as well as the doped states 1,2 . 
There are also several magnetic compounds which can be considered as coupled 
spin cluster systems in which the dominant exchange interactions occur within 
clusters of spins. The clusters are coupled through weaker exchange interactions. 

SA prominent example of such systems is that of molecular magnets 3 . Examples 
of spin clusters are dimers and four-spin plaquettes. The magnetic properties 
*^J ■ of spin clusters can be determined exactly if the size of the cluster is small. It 

£h ! is of significant interest to find out how the cluster properties are modified in 

the bulk. In this paper, we propose a class of two-chain spin ladder models 
which can be defined in terms of n-spin clusters and dimers. The specific values 
of n considered are n = 4 and 6 though a generalisation to higher n values 
is also possible. The models describe spins of magnitude 1/2 interacting via 

jJJ ■ modulated antiferromagnetic (AFM) exchange interactions. We show that in 

a wide parameter regime, the exact ground and low-lying excited states of the 
full ladder model can be described in terms of the eigenstates of the dimers 
and n-spin clusters, i.e., the clusters act as decoupled entities even in the bulk. 
This is also true when the exchange interactions coupling the clusters are of 
considerable strength. In the same parameter region, the magnetization/site of 
the ladder in the presence of an external magnetic field exhibits the phenomenon 



of magnetization plateaus. The condition for the appearance of a plateau is given 
by 4 

S u — m-u = integer (1) 

where S u and m u are the total spin and magnetization in unit period of the 
ground state. 

Spin-1/2 ladder models of various types have been extensively studied both 
in zero and finite magnetic fields 5-16 . In frustrated spin ladder models, extra 
diagonal exchange couplings (one or two) are present. Hakobyan et al. 10 have 
given an overview on the phase diagram of the general frustrated two-chain 
ladder model. The magnetization process of the general model is, however, yet 
to be undertaken. Brenig et al. 11 have introduced a dimerized and completely 
frustrated two-chain ladder model. The model, in which diagonal exchange in- 
teractions are present in every plaquette, is equivalent to a chain of edge sharing 
tetrahedra. The rung exchange interactions in the ladder model are of strength 
Ji. The nearest -neighbour intra-chain and diagonal exchange interactions are 
of equal strengths and in alternate plaquettes the strengths are J\ and J3 re- 
spectively. Examples of tetrahedral cluster compounds are tellurates of the type 
Ca^FeiO^Xi with X — CI, Br 12 . The class of models we propose describes 
two-chain spin-ladders with modulated exchange interactions. The spin ladders 
consist of tetrahedral spin clusters containing n spins separated by two-spin 
clusters, i.e., dimers (rungs). In section II, we introduce the models and dis- 
cuss the magnetization properties in the presence of an external magnetic field. 
Section III contains a summary and discussion of the major results obtained. 

The spin ladder model consists of blocks of tetrahedral clusters containing 
n spins separated by rungs of two spins. Fig. 1. shows the simplest such ladder 
model with n = 4. The tetrahedral clusters are represented by solid lines and 
the two-spin rungs (dimers) by dashed lines. Within a tetrahedral cluster, the 
diagonal exchange interactions are of strength J3 and the other exchange inter- 
actions are of strength J\ . The rung exchange interaction strengths are J' and 
a rung is coupled to a neighbouring tetrahedral cluster through exchange inter- 
actions (dashed lines) of strength J 2 . Periodic boundary condition is assumed 
to hold true. The spin Hamiltonian describing the ladder model is given by 

H= 2_^ [Jl(S \i. S u+i + S 2i- S 2J+1 + S lj. S 2i + S lj+l. S 2i+l) 

t=3j+l,j=0,l,- 

+ Jz{S li- S 2i+l + S 2i- S li+l)\ + J 2_^ £>li-&2i 

i=3j,j=0,l,- 

+ J2 2_^i (S u + S 2i + S u + 2 + S 2i+2)-(S u+i + S 2i+l) (2) 

i=3j+2,j=0,l,- 



H T + H R + H T r 
2 



The spin operator S u ( S ii ) is associated with the z-th site of the lower 
(upper) chain of the ladder, the site indices are sequential as shown in Fig. 
1. The Hamiltonians Ht and Hr describe the tetrahedral clusters and the 
rungs respectively whereas Htr contains the exchange couplings between the 
tetrahedral clusters and the rungs. The total spin of each rung is a conserved 
quantity due to the special structure of the Hamiltonian. 

We now determine the ground state of the ladder model. Using the method 
of 'divide and conquer' 17 , it is easy to show that for J2 < 4~' trie exact ground 
state has all the rungs and the tetrahedral clusters in their ground state spin 
configurations. A brief sketch of the proof is given in the following. The ground 
state of a rung is a singlet. The ground state of a tetrahedral spin cluster is 
a resonating valence bond (RVB) state and has total spin S=0. The ground 
state is \ipRVBi) for J3 < J\ and \4>rvb2) forJ3 > J\. The states \ipnvBi) and 
\1pRVB2) are linear combinations (plus and minus) of two valence bond (VB) 
states. In one VB state, the two VBs (singlets) are horizontal and in the other 
vertical. 
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Table I: The energy eigenvalues and eigenvectors of a tetrahedral cluster 
with exchange interactions of strength J\ (horizontal and vertical) and J3. The 
eigenstates \tpRVBi) and \1pRVB2) are the resonating valence bond states. 



At J3 = Ji, the ground state of a tetrahedral cluster is doubly degenerate. The 
two states have a pair of singlets (valence bonds) along either the horizontal or 
the vertical bonds. The Hamiltonian Htr (Eq. (2)), containing the exchange 
couplings between the rungs and the tetrahedral clusters, act on the rung singlets 
to give zero and thus has no contributions to the energy E\ of the eigenstate. 
E\ is the sum of the ground state energies of the tetrahedral clusters and the 
rungs. Let E g be the exact ground state energy of the total Hamiltonian H. 
Then E g is < E\. Let \tp g ) be the exact ground state wave function. Then, 
from variational theory, 



E g = j2 ^ a \ H i i^> + E Wfci H 'i i^> ^ E(^° + E 'io) 



(3) 



H = Y,{Hi + Hi) 



where Hi 's are the tetrahedral cluster Hamiltonians with the ground state 
energy Ei (Table I) and H^s are the six spin cluster Hamiltonians, each of which 
contains the rung exchange interaction Hamiltonian and the eight exchange 
couplings (four horizontal and four diagonal) which connect the rung to nearest- 
neighbour tetrahedral clusters. The ground state energy of H[ is E' io . For 
Ji < 4"' ^'io 1S tne ener §y °f a singlet across the rung. We can now write down 
the inequality, 

Y,{E io +E' io )<E g <E l (4) 

i 

E\ is, however, exactly equal to ^2n{Ei +E' io ) since it is the sum over the ground 
state energies of all the rungs and the tetrahedral clusters. Thus, E g — E\, i.e., 
the exact eigenstate is also the exact ground state of the full ladder model. The 
ground state has the novel structure of islands of four-spin RVB configurations 
in the tetrahedral clusters separated by singlet spin configurations along the 
rungs. The exact ground state energy is E g — N(Em — 3^-) where N is the total 
number of tetrahedral clusters as well as rungs in the ladder. E i0 = —2J\ + 4g- 
for J 3 < Ji and E i0 = -^ for J 3 > Ji (Table I). When J 3 = Ji, the exact 
ground state is highly degenerate. The number of such states is 2 N . 

We now want to check whether the exact ground state is still a product of 
the ground states of the rungs and the tetrahedral clusters when J 2 is made 
larger than ^ . For this, the total Hamiltonian H (Eq. (2)) is written as a sum 
over six-spin sub-Hamiltonians, hi's, i.e., H = X^^»- Each sub-hamiltonian 
describes a tetrahedral cluster coupled to a rung. The six-spin sub-Hamiltonian 
can be diagonalised exactly to obtain the ground state energy. Again, one uses 
the method of 'divide and conquer'. When the six-spin sub-Hamiltonians are 
added together to obtain the full Hamiltonian, the J\, J 3 , J' bonds are counted 
twice and the Ji bonds only once. One can identify the region of parameter 
space in which the exact ground state of the full ladder is of the product form. 
Fig. 2 shows the phase boundaries, in the parameter space of j 2 - and y- for 
different values of -j 2 - . In the parameter regime below each phase boundary, 
the exact ground state is a product over the ground states of the rungs and 
the tetrahedral clusters. One finds that in certain parameter regimes Ji can be 
larger than ^ and the exact ground state continues to be of the product form. 

We next include an external magnetic field term — /i^i=i &i m the Hamil- 
tonian H (Eq. (2)), where 6N is the total number of sites in the ladder. We first 
consider the case of a single tetrahedral cluster in the presence of a magnetic 
field. The magnetic field couples to the z-component of the total spin of the 
cluster, S^ ot , which is a conserved quantity. The ground state energy E g {S^ ot ) at 
h = for Sf 0t = 0,1 and 2 can be obtained from Table I. When the external field 
h ^ 0, the ground state in each S% ot subspace is E g {Sf ot , h) = E g {Sf ot ,Q) — hSt ot . 



The ground state magnetization curve can be easily obtained. Consider the case 
J3 < J\. The magnetization per site m is zero from h = upto a critical field 
h ci = J\- For h ci < h < h C2 — 2J l5 ra=\ and beyond h — h C2 , the saturation 
magnetization, m — i, is obtained. Thus there are three magnetization plateaus 
at m = 0, j and \. For the external field h = 0, we have already seen that there 
is an extended parameter regime in which the exact ground state of the full 
ladder is a product of the ground states of the rungs and the tetrahedral clus- 
ters. We now investigate whether the same holds true for a finite magnetic field. 
Again, one uses the method of 'divide and conquer' and the sub-Hamiltonian 
used is a six-spin cluster consisting of a tetrahedral cluster and a rung. For the 
full ladder, one can identify a region (region A) in parameter space in which for 
< h < h Cl , in is zero. At h Cl , there is a jump in the value of m to m = | and 
a plateau is obtained for h upto h C2 (Fig. 3). When h Cl < h < h C2 , the exact 
ground state has the tetrahedral clusters in their S z = 1 ground states and the 
rungs in singlet spin configurations. Since, the number of tetrahedral clusters is 
N and the total number of sites is 6N, the magnetization/site m in the ground 
state is g. The quantization condition in Eq. (1) is obeyed as unit period of 
the ground state contains six spins so that S u = 3 and the magnetization m u 
in the unit period is 1. At h C2 , there is a second jump in m from g to |. When 
h C2 < h < h C3 , the exact ground state has the tetrahedral clusters in their 
S z — 2 ground states and the rungs in singlet spin configurations. In this case, 
S u and m u in Eq. (1) are 3 and 2 respectively. At ft. = h C3 , there is a jump in m 
from I to the saturation magnetization \. For J3 < J\, h Cl , h C2 and h C3 have 
the values J\, 2J\ and J' + J2, (2Ji < (J' + J2)) respectively. There are other 
parameter regions (regions B and C) in the parameter space in which the full 
plateau structure in the m versus h plot, as shown in Fig. 3, is not obtained. 
Fig. 4 shows the phase diagram for the full ladder in a magnetic field in the -j- 

vs. -j 2 - parameter space and for j 2 - = 0.2. The region A exhibits the full plateau 
structure in m vs. h as shown in Fig. 3. In region B, the jump in m from to | 
occurs at h = h Cl (Fig. 3) but beyond h C2 , the ground state is no longer of the 
product form. In region C, the ground state loses its simple product structure 
beyond h = h Cl . Similar phase diagrams are obtained for higher values of -j^-and 
also for J3 > J\. One can generalise the ladder model shown in Fig. 1 by as- 
signing different coupling strengths J\, J a and J3 to the vertical, horizontal and 
diagonal couplings of the tetrahedral clusters. Again, results similar to the case 
3\ = J4 are obtained. With J4 = J3 = J2 and J\ = J' , the two-chain frustrated 
ladder model introduced by Bose and Gayen 18 is recovered. In a finite magnetic 
field h, the magnetization/site m vs. h has a simple plateau structure 6 . 

Another generalisation of the ladder model shown in Fig. 1 is to replace a 
tetrahedral cluster by a block of tetrahedral clusters. Fig. 5 shows an example 
in which the block contains two tetrahedral clusters. The six-spin blocks are sep- 
arated by two-spin rungs. Again, one can show that in an extended parameter 
regime, the ground state has the product form in both zero and finite magnetic 
fields. The exact ground state is the product of the ground states of the six-spin 
blocks and the rungs. The ground state of a six-spin block is a RVB state. An 



extra magnetization plateau exists for h C3 < h < h Ci in which the ground state 
has all the six-spin blocks in their S z = 3 ground state configurations and the 
rungs are in singlet spin configurations. At h = h C4 , m jumps to its full satu- 
ration value. Fig. 6 is the phase diagram similar to Fig. 4 for the full ladder 
with jp- = 0.2. In the 'divide and conquer' method, the full ladder Hamiltonian 
is a sum over eight-spin sub-Hamiltonians. Each sub-Hamiltonian describes the 
interactions in a block of spins consisting of two tetrahedral clusters and a single 
rung. In region A\, the full plateau structure in m vs. h is obtained. In regions 
B\, C\ and D\, the ground state no longer has the product form beyond the 
fields h C3 , h C2 and h Cl respectively. Similar phase diagrams are obtained for 
higher values of y-. One can generalise the ladder models shown in Figs. 1 and 
6 by making the blocks of tetrahedral clusters of bigger size (the total number of 

spins in a block may be 4, 6, 8, 10 etc.). Two-spin rungs separate the blocks 

of spins. In certain parameter regimes, the exact ground state is possibly the 
product of the exact ground states of the rungs and the blocks of tetrahedral 
clusters. A full study of such ladder models is yet to be undertaken. 

In this paper, we have described a class of two-chain ladder models consisting 
of blocks of tetrahedral clusters, containing n spins, separated by two-spin rungs. 
We have specifically considered two cases: n=4 and 6. We have shown that in 
an extended parameter regime, the ground state of the ladder is a product over 
the ground states of the rungs and the blocks of tetrahedral clusters. For n=4, 
we have shown that the exact ground state consists of RVB spin configurations 
in the tetrahedral clusters and the rungs are in singlet spin configurations. For 
J3 = Ji, the ground state is highly degenerate. When J3 is < J\ (> Ji), the 
tetrahedral cluster is in the RVB state \tpRVBi) (\1pRVB2)) and the exact ground 
state of the full ladder model is non-degenerate. A notable feature of the ladder 
model is the presence of singlet excitations in the triplet spin gap in certain 
parameter regimes. As already pointed out in earlier references 11 ' 12 , the singlet 
energy level \1pRVB2) of a tetrahedral cluster crosses the triplet energy level at 
J3 = 4r (Table I). Thus for 4p < J3 < J\, the singlet excitation described by 
\1pRVB2) falls in the triplet gap. Similarly, for J 3 > J\, \1pRVB2) is the ground 
state and for J\ < J3 < 2Ji, the singlet excitation corresponding to \?Prvbi) 
falls in the triplet gap. These features carry over to the case of the full ladder 
model in the parameter region in which the exact ground state can be written 
in a product form. The existence of singlet excitations in the triplet spin gap 
is a characteristic feature of some other AFM spin systems which include the 
S = 2 Heisenberg antiferromagnet (HAFM) on the kagome lattice 19 , the S = 5 
HAFM on the pyrochlore lattice 20 and some S = \ AFM spin models on the 
^-depleted square lattice 21 ' 22 . 

The model shown in Fig. 1 can be generalised to bigger blocks of tetrahe- 
dral clusters. (Fig. 5 shows blocks of two tetrahedral clusters). Instead of the 
tetrahedral cluster shown in Fig. 1, one can also consider a generalised tetrahe- 
dral cluster with the horizontal, vertical and diagonal exchange interactions of 
different strengths. Again, in an extended parameter regime, the exact ground 
state is found to be of the product form. The ground states have the interesting 
structure of islands of RVB spin configurations separated by singlet spin con- 



figurations along the rungs. This type of exact ground state is not known for 
other spin models including ladders with modulated exchange interactions. 

The ladder models have also been studied in an external magnetic field h. In 
the parameter regime in which the ground states in the different magnetization 
subspaces are of the product form, the magnetisation/site m as a function of h 
exhibits plateaus (Fig. 3). The quantization condition in Eq. (1) is obeyed at 
each plateau. Figs. 4 and 6 show the phase diagrams for the ladder models of 
the types shown in Figs. 1 and 5. Both the phase diagrams show that there 
are extended regions in parameter space in which the ground states in differ- 
ent magnetization subspaces are of the product form. Kolezhuk 23 has studied 
magnetization plateaus in a spin system consisting of strongly coupled dimers 
which are again weakly coupled in a planar arrangement of zigzag interactions. 
In our ladder models, we have two different type of clusters: dimers (two-spin 
rungs) and tetrahedral clusters. Further studies are needed to obtain the phase 
diagrams of the ladder models in the full parameter space. 

The authors thank S. Ramasesha and K. Tandon for letting them use their 
Heisenberg Calculator (exact diagonalisation program) for quantum spin sys- 
tems. E. Chattopadhyay is supported by the Council of Scientific and Industrial 
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Figure Captions 

Fig. 1. Two-chain ladder model consisting of tetrahedral clusters (solid lines) 

coupled to two-spin rungs (dashed lines). The exchange interaction strengths 

are as shown in the Figure. 

Fig. 2. Phase diagram of the ladder model (Fig.l) in the parameter space of 

■j 2 - and -j- . The parameter space below a solid line corresponds to the phase in 

which the exact ground state is a product over the ground states of the rungs 

and the tetrahedral clusters. 

Fig. 3. Plot of magnetization/site m versus external magnetic field h for the 

two-chain ladder model shown in Fig. 1. The plot is obtained in the parameter 

region in which the exact ground states in different Sf ot subspaces have the 

product form. Two non-trivial magnetization plateaus occur at m = \ and 

m=\ . 

Fig. 4. Phase diagram of the ladder model (Fig. 1) in a finite magnetic field 

and in the parameter space of j- and ^f- with j 2 - = 0.2. The regions A, B and 

C are explained in the text. 

Fig. 5. A two-chain spin ladder which consists of blocks of two tetrahedral 

clusters coupled to two-spin rungs (dashed lines). The exchange interaction 

strengths are as shown in the Figure. 

Fig. 6. Phase diagram of the ladder model (Fig. 5) in a finite magnetic field 

and in the parameter space of jj- and j^- with j^ = 0.2. The regions A\, B\ , 

C\ and D\ are explained in the text. 
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